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Abstract. We propose a new theory for the formation of rRi ring structures, i.e. for ring structures with both an inner and an 
outer ring, the latter having the form of "8". We propose that these rings are formed by material from the stable and unstable 
invariant manifolds associated with the Lyapunov orbits around the equilibrium points of a barred galaxy. We discuss the 
shape and velocity structure of the rings thus formed and argue that they are in agreement with the observed properties of rRi 
structures. 
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1. Introduction 

Barred galaxies often show spectacular rings whose different 
types have been classified by Buta ( 1986ai as follows: Nuclear 
rings (which are not discussed here) which surround the nu- 
cleus and that are much smaller in size than the bar. At larger 
radii there are inner rings, denoted in Buta's ( 1986ai classifica- 
tion by a "r", which suiTound the bar and which have the same 
size and orientation as the bar. And there are outer rings, de- 
noted by R, which are bigger than the bar. These "pure" rings 
are defined to be distinct and closed, but one can often find 
unclosed or partial ring patterns of spiral character, and these 
are referred to as "pseudorings" and denoted R'. A particular 
class of outer rings called Ri or Ri', for pseudorings, has two 
main arms forming an eight-shaped ring or pseudoring, with its 
major axis perpendicular to the bar. NGC 1326 is a well stud- 
ied example of an {Ri)S B{r)0/a galaxy (Buta , 1995) . Its bar is 
surrounded by an inner ring which is almost exactly aligned 
with the bar and has roughly the same diameter. Its outer ring 
is clearly Ri and is elongated perpendicular to the bar (Buta et 
al. rT998l see Fig.©. 

Schwarz ( flMTl TM^ [T985I showed that ring-like struc- 
tures can arise around the Lindblad resonances due to a bar-like 
perturbation of the galaxy potential. The gas will be forced to 
reaiTange its distribution and generate a spiral. Near the outer 
Lindblad resonances the crossing of perturbed trajectories will 
develop a ring-like pattern. In these regions, gas clouds will 
collide and will form spiral shock fronts which will slowly 
change as a result of torques exerted by the bar and evolve into 
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a ring structure which, after star formation, will be populated 
by stars in near-resonant periodic orbits. 
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II 9651 1 argued that orbits in the gravitational poten- 
tial of a bar play an important role in the formation of arms. He 
noted that orbits departing from the vicinity of the equilibrium 
points located at the ends of the bar describe loci with the shape 
of spiral arms and can be responsible for transport of stars 
from within to outside corotation, and viceversa. Unfortunately, 
he did not set his work in a rigorous theoretical context, so 
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that it remained purely phenomenological. He also investigated 
whether orbits can be responsible for ring-like structures, but, 
in this case, he did not consider orbits departing from the ends 
of the bar as he previously did when accounting for the spi- 
ral arms. He concluded that rings would require high energy 
orbits, while we will see in this paper that mainly low energy 
orbits can constitute the rings. 

Here we propose a new dynamical model, applicable to the 
particular case of rRi class of ringed galaxies. We expect that 
more detailed modelling will extend this new model to the rest 
of the ringed galaxy classes. The model is based on the orbital 
motion of stars in the vicinity of equilibrium points in the ro- 
tating bar potential and does not rely on additional star forma- 
tion. The tools we use in our model are well known in celes- 
tial mechanics but have not been much used so far in galactic 
dynamics. We therefore start, in section 2, by introducing the 
equations of motion and the equilibrium points. In section 3 we 
describe the dynamics around the Lagrangian point Li and we 
introduce the invariant manifolds. We discuss the linear case in 
subsection 3.1 and the general case in subsection 3.2. In section 
4 we describe the role of the invariant manifolds in the transport 
of stars and their properties in the framework of ringed barred 
galaxies. In section 5 we describe briefly how invariant mani- 
folds are computed numerically and we apply them to a specific 
model. Finally, section 6 discusses the properties of the spiral 
arms thus produced and summarises our results. 

2. Equations of motion and equilibrium points 

We model the potential of the galaxy as the superposition of 
two components, one axisymmetric and the other barlike. The 
latter rotates clockwise at angular velocity £lp - QpCz, where 
Qp > is the pattern speed considered here to be constant ' . 
The equations of motion in the frame rotating with £lp are: 

f = -VO - 2(ap X f ) - £lp X (Hp X r), (1) 

where the terms -2£lp x r and -Hp x (Qp x r) represent the 
Coriolis and the centrifugal forces, respectively, and r is the 
position vector 

Following Binney & Tremaine (fl987% we take the dot 
product of equation Q with f , and rearranging the resulting 
equation, we obtain 

dEj 
dt 



0, 



where 



E, = 



-H<D- 



1 



fip xr I 



equation Q becomes 

f = -VOeff - 2(Qp X f), 

and the Jacobi constant is 



(2) 



Ej^-\r\^+ (D, 



eif, 



SO that it can be considered as the "energy" in the rotating 
frame. 

The surface Ogff - Ej \?, called the zero velocity surface 
and its cut with the z = plane is the zero velocity curve. All 
regions in which Oeff > £7 are forbidden to a star and we call 
them forbidden regions. Fig.O> shows an example of zero ve- 
locity curves and the regions delimited by them, namely the 
exterior, interior and forbidden regions, for the potential intro- 
duced in section 5. 

<l>eff has five equilibrium points, named L\ to L5, located 
in the xy plane, at which - = = 0. Due 

to their similarity to the corresponding points in the restricted 
three body problem, they are often called Lagrangian points. 
L\ and L2 lie on the x-axis and are symmetric with respect to 
the origin. L3 lies on the origin of coordinates. Finally, L4 and 
L5 lie on the y-axis and are also symmetric with respect to the 
origin (see Fig|2t, again for the potential introduced in section 
5). 

We can check the stability of the Lagrangian points by con- 
sidering the motion in their immediate neighbourhood^. If we 
expand (Jeff around one of these points and retain only first or- 
der terms, the equations of motion become 



( X = 2Q.py - <^x^x 
\ y = -2Qpi - OvvJ 

where we have defined 



(3) 



x = x-X]^; y = y-yi^; z = Z-Zi^, 



eif 



&2 



Ej known as the Jacobi integral or Jacobi constant. Notice 










that this is the sum of ^f^ -1- which is the energy in a nonro- 




= fiixu. 


. . , xe) 


— X4 


tating frame, and of the quantity 


X2 


= fiixu- 


. . , xe) 


= X5 


-5 1 Qp X r (x^ + y^), which can be fliought of as the 


x-i 


= Mxu- 


. . , xe) 


- Xe 


"potential energy" to which the centrifugal "force" gives rise. 


X4 


= f4(xi,. 


. . , xe) 


= 2Qp X5 - (b^;,Xl 


So if we define an effective potential 


X5 


= Mxi,. 


. . , Xe) 


- -2Q.p X4 - <S>yyX2 


,1 „i , 1 9, 




= Mxi,. 


. . , Xe) 


= -^7zX3- 



and Xl, and zl are the coordinates of the Lagrangian point 
Lj. Note that, for any barlike potential whose principal axes lie 

along the coordinate axes, ^ ''^oy ) = by symmetry. Setting 

Xl — X, X2 — y, x-i — z, X4 = X, X5 = j and xe = z, 
equations ^ are written as a system of first order diff'erential 
equations. 



(4) 



Oetf = D ■ 



Bold letters denote vector notation. 



- We can refer to Pfenniger H990'\, where he studied the stability 
character of the Lagrangian points in different stellar bars. 
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Fig. 2. a) Outline of the bar and the position of the five equilibrium points (marked with a star), b) The two Lyapunov orbits (thin 
solid lines) and the zero velocity curves (dot-dashed lines) delimiting the forbidden, the interior and exterior regions. Both are 
given for the same energy value. 
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3. Dynamics around L\. The invariant manifolds 

3A. Linear case 

Let us now focus on the dynamics around Li (L2 is completely 
symmetrical). The differential matrix associated to system (0) 
around Li is 



Dfx(Li) = 



We obtain the stability character of Li by studying the eigen- 
values of this matrix. It has six eigenvalues: A, -A, toi, -a>i, vi 
and -vi, where A, co and v are positive real numbers, i.e. Li is 
a linearly unstable point. The corresponding eigenvectors (ei- 
ther real or complex) have zero X3, and components in the 
case of +A and +a)i, and zero xi , X2, X4, and xs components in 
the case of ±vi. Since the purely imaginary eigenvalues denote 
oscillation and the real eigenvalues are associated to a saddle 
behaviour, i.e. exponential behaviour with opposite exponents, 
the linearised flow around Li in the rotating frame of coordi- 
nates is characterised by a superposition of an harmonic mo- 
tion in the xy plane (equatorial plane), a saddle behaviour in 
this plane, and an oscillation in the z-direction. 

Because of this unstable character, the equilibrium points 
Li and L2 set the limits of the stability region around the stable 
point L3, and thus, they set an upper limit to the extension of the 
bar. This, however, is only an upper limit and in most cases the 
bar is much shorter than that limit ( Athanassoula . 1 992m Patsis, 
Skokos & Athanassoula 2003). 

The central stable point, L3, is suiTounded by the classic 
xi family of periodic orbits which is responsible for maintain- 
ing the bar structure, while the stable points L4 and L5 are sur- 



rounded by families of periodic banana orbits (Contopoulos & 
Papavannopoulos f 1 9801 Athanassoula et al. 1983 ; Contopoulos 
11981, Skokos, Patsis & Athanassoula 2002). All these orbits 
have been well studied for many models (Contopoulos 120021 
and references therein) and we will not discuss them here any 
further. 

As already mentioned, the general linear motion around Li 
is obtained by the addition of an hyperbolic exponential part to 
the in-plane and out-of-plane oscillations mentioned. We note 
that this exponential part has both stable and unstable compo- 
nents with exponents of opposite sign. Following again Binney 
& Tremaine J1987> we write 



(5) 



x(t) = Xie-^' + Xie'-^' + X3 cos(dL»f + 0), 
y{t) = X^e'^' + X^e^'^' + Xf, sin(wf + if>), 

z(t) - X7C0S(vf -I- (f/). 



Here Xj,i - 1, . . . , 7 and 0, t// are values representing ampli- 
tudes and phases. Substituting these equations into the differ- 
ential equations (|3|l, we find that X, are related by 



X4 
X5 



2n,,A 



■Xi = 



2£V V 



2n„A 



■X2 = 



0„-a;- y 



<ti„+A 
20. 



X2 



(6) 



-^3 



Note that A 1 depends 



We define Ai = %r4- and Ar - 
only on A and A2 depends only on oj. Moreover Xi, - A\Xi, 
Xi - -A1X2, and X(, = A2X3, so that equation (|5} becomes 



x(t) = Xie-^' + X2e''^' + X3 cos(wf + (p), 

y(t) = AiXy - AxXte-^' + A2X3 sin(wf -1- 0), 

z(i) - Xj cos(vr -I- f ). 



(7) 
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In the sequel we will restrict ourselves to the motion in the 
equatorial plane (i.e. z - 0, or X-j = 0). This restriction is not 
critical in the dynamics we want to study, since the z com- 
ponent essentially only adds a vertical oscillation to the pla- 
nar motion. For a study of the vertical orbital structure around 
the Lagrangian points in barred galaxies, see Olle & Pfenniger 
iT998l . 

Using Q, a given state (x, y, i, j) at f = is characterised 
by a choice of (Xi,X2,X3,0i), modulus In in the phase (p. 
When Xi = X2 - the initial condition is of the form 

(X3 COS (p, A2X3 sin (p, -XiO) sin 0, cos (p) 

for selected values of X3 and (p. When time evolves, we obtain 
from this initial condition the periodic motion, 

Xo(f) = (x,y,x,y) = 

(X3 cos(dLtf + tp), A2X3 sin(dLtf + tp), 
-XjCl) sin(wf -I- (p), A2XT,a> cos{a>t + (p)), 

of period t to which we will refer as a linear Lyapunov periodic 
orbit. 

Consider now any small deviation 6 from the periodic orbit 
xo(f), x(f) - X()(f) + 6. Inserting this into the equations of mo- 
tion 0, and linearising them with respect to 6, we obtain the 
variational equations 



as defined in Contopoulos J2002> . where 

(dff 



(8) 



A(f) = {6ik) 



dxi 



( ££l 



3.1-6 



is the variational matrix. The variational equations are linear 
equations in 5, (/ = 1, . . . , 6) with periodic coefficients of pe- 
riod T. These equations are used to study the stability character 
of a periodic orbit. If we integrate the variational equations (|8|i 
until time t with initial conditions 6\k - (1, 0, . . . , 0), . . . , = 
(0, . . . , 0, 1), we obtain the monodromy matrix. The eigenval- 
ues and eigenvectors of the monodromy matrix give us infor- 
mation on the stability character of a periodic orbit. 

Returning to equations (0, let us consider a similar initial 
condition but with X\ - ^ and X^ + 0. According to the 
exponential term proportional to X^ vanishes when time tends 
to infinity and the trajectory tends to the linear Lyapunov. All 
these type of orbits form what is called the stable manifold of 
the linear Lyapunov. In the same way, if the initial condition is 
chosen with X\ 0, X2 = the exponential term proportional 
to X\ tends to zero when time tends to minus infinity, and all 
such orbits form what is called the unstable manifold of the lin- 
ear Lyapunov. Roughly speaking, orbits in the stable/unstable 
manifold are asymptotic orbits, which tend to/depart from the 
linear Lyapunov orbit. 



Initial 

Condition on 
Unstable Invariant 
Manifold (f/),0 




Lyapunov orbit 



Unstable 
Invariant 
Manifold 



Fig. 3. Schematic view of an unstable invariant manifold. 
3.2. General case 

All the definitions given in section 3.1 for the linear case are 
easy to extend to the general case when the full equations of 
motion are considered. From Li emanates a family of planar 
periodic orbits known as Lyapunov orbits (Lvapunov I1949> 
which locally can be parametrised by the energy. The full set 
of Lyapunov orbits form a family of periodic orbits, which we 
will denote by F. 

Let us denote by ^^(f , X) the orbit that at f = has the state 
X = (x,y,x,y). For a given y e F, i.e. for a given Lyapunov 
orbit, we define the stable manifold of y as 



w: 



/lim||»l'(f,^)-rll=0 

/ — >oo 



where the double bars denote Euclidean distance. Thus, simply 
speaking, the stable manifold is the set of orbits which tend to 
the Lyapunov orbit as time tends to infinity. In the same way, 
the unstable manifold of y is defined as 



/ lim ||>I'(f,^)-rll=0 



Simply speaking again, the unstable manifold is the set of or- 
bits which tend to the Lyapunov orbit as time tends to minus 
infinity, or, equivalently, the set of orbits departing from the 
Lyapunov. The orbits of and W" have the same energy as 
the Lyapunov orbit y, and so they belong to the same energetic 
three dimensional manifold where y is contained. Moreover, 
Wy and W" are two dimensional tubes, which, similarly to the 
former linear example, can be parametrised by the angle <p and 
the time f (see Fig.l^and Masdemont ( 2005 )). We also note that 
both Wy and Wy have two branches meeting at the Lyapunov 
orbit in a way similar to a saddle point (see Fig. Ell. 



4. The role of invariant manifolds 

Invariant manifolds of Lyapunov orbits play a crucial role in 
the transport of material between different parts of the config- 
uration space. Lyapunov orbits are located near the ends of the 
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Fig. 4. Invariant manifolds. In the centre of the plot, and 
in white solid line, the Lyapunov orbit around Li. The two 
branches of the unstable invariant manifold (red or dotted 
lines), and the two branches of the stable invariant manifold 
(green or dotted lines). In grey, the forbidden region surrounded 
by the zero velocity curves. 



bar of the galaxy, between the two banana like zero velocity 
curves which surround the forbidden regions of motion for the 
considered energy, and, loosely speaking, can be considered as 
gates between the interior and exterior region they delimit (see 
FigO>). Let us consider for instance the stable manifold, W*, 
inside the interior region and integrated backwards in time till it 
crosses the plane, S , defined by x = 0. The plot of this intersec- 
tion in the yy plane is a closed curve j (Fig.|5}. Each point 
in the plane yy of S corresponds to a given trajectory, since 
.X = by the definition of S and x can be obtained from the 
condition that the energy of the state (x, y, i, y) be the selected 
one (the sign of i is determined by the sense of crossing). j 
is a closed curve that splits the yy plane in S in three differ- 
ent regions: the curve itself, the points exterior to the curve 
and the points interior to the curve. By definition, the points on 
the curve , belong to W:' and are therefore orbits that tend 
asymptotically to the Lyapunov orbit. The points outside j 
are states whose trajectories remain inside the interior region 
of the galaxy delimited by the zero velocity curves, while the 
points inside j coiTespond to orbits that transit from the in- 
terior region to the exterior one. These last orbits, the transit 
orbits, are confined inside the tube W" and, as we will argue 
in the following sections, are the orbits which form part of the 
rings of the galaxy for the considered energy value. In this way, 
the manifolds of the Lyapunov orbits drive the motion of the 
stars from the interior to the exterior regions. For more details 
of this mechanism, although in another context, see Gomez et 
al. (2004) and references therein. Since these invariant mani- 
folds are not limited to the vicinity of the unstable points, but 
extend well beyond it, they can be responsible for global struc- 
tures and we will argue in this paper that they, together with the 
orbits driven by them, could be responsible for the ring struc- 
tures in barred galaxies. 




-4.6 -4,5 -4,4 -4,3 -4,2 

y 




Fig. 5. Transport of material. Top panel: the curve j in the 
yy plane. Bottom panel: schematic view of the dynamics in the 
region around the Li Lagrangian point. In grey, the region de- 
limited by the zero velocity curves. The dotted line gives a non- 
transit orbit, which is confined to the interior region and whose 
intersection with the yy plane is located outside the W'^ j curve; 
the solid line, a transit orbit, which in the yy plane is located 
inside the j curve; and the dot-dot-dash line, an asymptotic 
orbit of the stable invariant manifold. In the inlay, a schematic 
view of the curve j in the yy plane with the location of the 
three orbits. 

5. Application to a ringed barred galaxy model 

In this section we will explicitly calculate the invariant man- 
ifolds in a barred galaxy model. The reader who has only 
skimmed the previous sections should keep in mind that the 
invariant manifolds are just ensembles of orbits linked to the 
Li and L2 Lyapunov orbits. 

Since invariant manifolds can only be calculated numer- 
ically, we first adopt a simple, yet realistic, barred galaxy 
model (Pfenniger" 1984V Our model bar consists of an axisym- 
metric component, modelled by a Miyamoto-Nagai potential 
(Miyamoto & Nagai[T975} 



■yJx^+y^ + (A+ Vb2T?)2 
and a Ferrers bar (Ferrers fTSTTt 

_ f pc(\ - m^y m < 1 

^ ~ I OT > 1, 

where = x^/a^ + f/b^ + z^/c^ and p, = is the 

central density. We take A = 3, B - 1, « = 2, a - 6, b - 1.5, 
c - 0.6, GMd - 0.9 and GMb =0.1. The pattern speed is taken 
such as to place corotation at the end of the bar. The length 
unit is the kpc, the total mass G(Md + Mb) is set to be equal to 
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-HOB -SOB 



X I 

Fig. 6. Unstable invariant manifolds for four different times 
given in the lower right corner of each panel. T is the period 
of the bar rotation. 

1, and the time unit is 2 x 10^ yr. With these units, the value 
of the potential energy at L3, Li(L2), and L^iL^) is -0.31503, 
-0.19789 and -0.19456, respectively. As done already in the 
previous sections, we limit ourselves to the z = plane, since 
the instability we are interested in is contained to this plane 
and the z component only adds vertical oscillations, which are 
unimportant in this context. 

In this model, we computed the invariant manifolds numer- 
ically using an approach similar to the one of Gomez et al. 
<1993> . i.e. we use a linear approach to obtain the initial con- 
ditions of the orbits that constitute the invariant manifolds^. 
As previously mentioned, the linear motion around L\ and L2 
consists of an hyperbolic exponential part in the in-plane and 
of an out-of-plane oscillation. The exponential part has both 
stable and unstable components, which correspond to the sta- 
ble and unstable eigenvectors of the monodromy matrix of the 
Lyapunov orbits around L\ and L2. Therefore, we obtain the 
initial conditions for the stable and unstable invariant mani- 
folds shifting positions and velocities of the Lyapunov orbit by 
a small amount (10^^) in the direction given by the stable and 
unstable eigenvectors, respectively. The global extension of the 
manifold is then obtained by integrating numerically with a 
Runge-Kutta-Felhberg of orders 7-8. 

Figure |6l shows the evolution of the length of the unstable 
invariant manifolds with time. Here, T is the rotation period of 
the bar. Since the invariant manifolds consist of a set of orbits, 
their time evolution is that of an orbit. So Fig.|6lshows the or- 
bits composing the invariant manifolds integrated up to 3,4,5 
and 6 bar rotations, respectively, in order to check their evolu- 
tion with time. Note that they leave the vicinity of the Lyapunov 
orbits, and stay close to the zero velocity curves. When they get 

' Higher order approximations could be obtained using similar 
techniques as in Masdemont ( 2005 (, but with much more difficulty 
and regularity problems. 



to the X - axis, they bend and head towards the opposite side 
of the bar, completing the ring structure. Thus our calculations 
show that the apocenter of the orbits is near the x - Q axis. 
This is indeed where we would intuitively place them since the 
periodic orbits outside corotation are elongated perpendicular 
to the bar and since the zero velocity curves also reach their 
maximum distance from the center at jc = 0. Also note that, 
although the Lyapunov orbit is unstable, the orbits in the man- 
ifold stay in its close vicinity for at least three bar rotations. 
By four bar rotation periods arm stubs are formed, while at five 
the arms have a winding of about 37r/4. By six bar rotation pe- 
riods the manifold has reached the opposite end of the bar In 
other words, the growth is by no means linear, since growth 
is slow in the beginning and faster as time advances. Indeed, 
it takes more than three orbital periods to leave the vicinity of 
the Lyapunov orbit, the first half of the winding is developed in 
about 4.5r and the full winding in barely at 6T . After bT, the 
orbits on the manifold are in the vicinity of the Lyapunov orbit. 
These orbits do not leave this region in another new direction, 
but they follow the direction given by the invariant manifolds 
already formed. 

Relevant information for a value of the Jacobi constant 
{Ej - -0.1977) close to the one of the L\{L2) equilibrium 
point, is given in Fig. In all four panels, the bar is out- 
lined by a black dot-dashed line, while the positions of the five 
Lagrangian points are marked with a star, and the dark grey 
lines correspond to the zero velocity curves of this particular 
value of Ej. This means that any orbit with this energy start- 
ing outside these curves cannot enter within them. We have 
also plotted the two plane unstable Lyapunov periodic orbits 
around the unstable Lagrangian points L\ and Li (black solid 
lines). In Fig.^t we show the unstable invariant manifolds for 
this value of Ej. Each of these manifolds is composed of two 
branches, an interior branch, lying in the interior region, and 
an exterior branch, lying in the exterior region. Each of these 
branches can be thought of as an ensemble of orbits moving 
away from the Lyapunov orbit. Both branches lie near the zero 
velocity curves and, as shown in Fig. |6l their length increases 
with time until they reach the opposite side of the bar from 
which they emanated. The interior branch, when complete, out- 
lines well the loci of the inner rings in barred galaxies. The ex- 
terior branch, when complete, has a shape similar to that of the 
R\ outer rings. From a dynamical point of view, in the phase 
space these branches are seen like tubes that drive the dynam- 
ics. In Fig. we plot the stable invariant manifolds, for the 
same value of E j. Again there are two branches, an interior and 
an exterior one. Note also that the space loci of the stable and 
unstable manifolds is almost identical. There is, however, an 
important difference in that, for the stable manifold, the orbits 
filling these tubes will go towards the Lyapunov orbits, while 
for the unstable manifold they will go away from it. In Fig. 

we show four sets of orbits starting from the vicinity of the 
Lyapunov orbits and following the four tubes that constitute the 
two branches of the unstable manifolds. This figure illustrates 
how the invariant manifolds drive the dynamics, since we can 
see that initial conditions in the vicinity of the Lyapunov orbit 
will follow, due to its unstable character, a trajectory close to 
the unstable invariant manifolds. In Fig. 01, we represent the 
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Fig. 7. Invariant manifolds and perturbations in the plane 
barred potential described in section 5 and for Ej = -0.1977. 
a) Invariant unstable manifolds coming from Lyapunov orbits 
around the Lagrangian points Li and L2. b) Invariant stable 
manifolds coming from the same Lyapunov orbits, c) Orbits 
starting from initial conditions near the Lyapunov orbit, d) 
Orbits starting from initial conditions inside the j curve of 
the interior branch of the stable invariant manifold. 

trajectories with initial conditions inside the j curve of the 
interior branch of the stable invariant manifold, as explained 
in the previous section. This set of orbits follows the stable 
branch they emanate from, approaching the Lyapunov orbit, 
and it leaves the bar region following the exterior branch of the 
unstable invariant manifold. 

6. Discussion 

The time during which the orbit stays around the Lyapunov or- 
bit before outlining the outer or the inner ring depends on the 
Lyapunov exponent (Lyapunov 1949 1 of the Lyapunov orbit, 
and is found to be an increasing function of Ej. Thus, orbits 
initially near the Lyapunov orbit with lower values of Ej stay 
less around the unstable point. This time increases gradually 
with the value of Ej, until the value at which the Lyapunov or- 
bit becomes stable after which all orbits starting in the vicinity 
of the Lyapunov orbit stay around it. Thus the orbits outlining 
the rings are mainly low energy orbits. Repeating the compu- 
tations in section 5 for different values of the Jacobi constant, 
Ej = -0.1973 and Ej = -0.1960, we find that the locus of 
the invariant manifolds (and therefore of the orbits associated 
to it) is roughly independent of their value of Ej. This is il- 
lustrated in Fig. Isji, for these two different values of Ej. As 
the energy increases, the size of the Lyapunov orbit also in- 
creases, so that the outline of the invariant manifold becomes 
thicker. However, as we consider more orbits and more energy 
levels, we find that the density of the central part of the out- 
lined area increases considerably, so that in practice the thick- 




Fig. 9. Velocity field along the invariant manifolds, a) In the 
rotating frame, b) In the nonrotating frame. 

ness of the ring will be considerably smaller than that of the 
higher energy manifolds. This is illustrated in Fig. [8^, where 
we plot the density profile on a cut across the ring. To obtain 
this figure we calculated the unstable invariant manifolds and 
the trajectories inside them for all energy levels at which the 
Lyapunov orbit is unstable, i.e. all energies for which the mech- 
anism we propose can be applied. This covers the range of val- 
ues from Ej = -0.19789, corresponding to the energy of the 
unstable equilibrium point, to the value Ej = -0.1674, where 
the Lyapunov family becomes stable. The contribution of each 
energy is weighted by a distribution function which is here sim- 
ply assumed to be an exponentially decreasing function of the 
energy, i.e. of the form exp (-^), with cr - 3Qkms~^ as a ve- 
locity dispersion characteristic of disc stars in the solar neigh- 
bourhood (Binney & Merryfield 1998). The exact shape of the 
distribution function is of little importance, but it has to be a de- 
creasing function of the energies in the rotating frame of refer- 
ence. The profile shown in Fig.|S^ is very similar to those found 
for the old and intermediate age stellar population in rings and 
spirals (Schweizer 1976). 

These ring structures will corotate with the bar. During 
bar evolution, however, the bar pattern speed will decrease 
with time due to an angular momentum exchange (Tremaine 
& Weinberg [T^ Weinberg 1985 Little & Carlberg 19^ 
r99Tb; Hernquist & Weinberg 1992; Athanassoula [19961 
Debattista & Sellwod I2OOOI Athanassoula .2003,. O'Neill & 
Dubinsky 2003; Valenzuela & Klvpin l2003t . This means that 
Li , L2, L4 and L5 will move outwards and thus the reservoir of 
fresh material for the rings would be continuously replenished. 
This may also be linked to the plumes that surround the ring 
structures of these galaxies (Buta 1984 1. 

The model we propose here forms ring structures which are 
not necessarily associated with the outer Lindblad resonance. 
In fact, in the model shown in Fig. the maximum radius of 
the inner branches of the invariants reaches values of about 4.25 
radial units, while corotation radius is placed at the end of the 
bar at 6 radial units. The outer branches extend up to 7.5 ra- 
dial units, while the Outer Lindblad Resonance of this model is 
placed at about 8.7 radial units. Note that the ratio of the maxi- 
mum radius of the outer and inner branches give a value of 1 .76 
which is in agreement with the values of this ratio measured for 
ringed galaxies (Athanassoula et al. 1 19821 Buta[^95 1. 

In the rotating frame, the velocities along the invariants 
reach a maximum at the point of maximum radius and then de- 
crease toward the region of the Lyapunov orbits (see Fig.|9^). 
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Fig. 8. a) Density profile on a cut across the ring, b) Two unstable invariant manifolds for different values of Ej. Note how similar 
the regions they delineate are. 



In the non-rotating frame the velocities along the invariants are 
an order of magnitude higher in the outer branches than in the 
inner ones (see Fig. |5J)). A clear prediction from our model 
is that the perturbations on the velocity fields produced by the 
rings should be higher around the outer ring than around the 
inner ring. There are little data available on the velocity fiels 
along the rings of Ri galaxies. The best data correspond to the 
galaxy NGC 1433 (Buta 1986b i, which is of type R\. For this 
case, we can tentatively say that there is a general agreement 
with the oscillations in the velocities measured in the rotating 
frame for the inner ring. 

Although these rings are not density waves, they do not 
have the shortcoming of material arms, since they do not wind 
up with time. A more appropriate name would be flux rings, 
since they are outlined by the trajectories of particles. The 
material of such rings would create a potential well. Other 
stars and gas in the galaxy would feel this potential and, while 
traversing the ring, they would stay longer at the potential min- 
ima, thus adding a density wave component to the ring. Thus, 
although we have not made any self consistent simulations, we 
can speculate that the flux rings and the density wave rings 
would coincide in the galaxy. 

In this paper we presented a new theory on the origin of 
rRi ring structures based on orbital dynamics. We introduced 
the invariant manifolds of a periodic orbit, used so far in ce- 
lestial mechanics. We explained their role in the transport of 
particles from the interior region to the exterior, and viceversa, 
and applied it to a realistic barred galaxy model. Finally, we 
compared and discussed our results and the characteristics of 
the rings obtained with observational data. We can conclude 
that the rRi ring structure can be interpreted as a bundle com- 
posed of all the invariant manifolds for all the possible energies 
as well as the orbits driven by them. 
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